Scaling laws for the response of nonlinear elastic media with implications for cell mechanics 
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We show how strain stiffening affects the elastic response to internal forces, caused either by material defects 
and inhomogeneities or by active forces that molecular motors generate in living cells. For a spherical force 
dipole in a material with a strongly nonlinear strain energy density, strains change sign with distance, indicating 
that even around a contractile inclusion or molecular motor there is radial compression; it is only at long distance 
that one recovers the linear response in which the medium is radially stretched. Scaling laws with irrational 
exponents relate the far-field renormalized strain to the near-field strain applied by the inclusion or active force. 
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The response of elastic media to internal forces is an impor- 
tant factor that governs the physics of defects and inclusions 
in solids and composite materials fjj]. Theory has shown that 
internal forces may be modeled as force dipoles that gener- 
ate strains and stresses that can lead to interesting collective 
effects [2]. Recently, similar ideas were applied to the me- 
chanical response of biological cells J3] where the resulting 
physical phenomena |5|,(5l] are also associated with cell func- 
tion H. 

The cytoskeleton of living cells contains molecular mo- 
tors that consume ATP and produce nonequilibrium forces by 
which cells attach to and pull on their surroundings. This me- 
chanical interaction relates to many aspects of cellular func- 
tion, from cell spreading and proliferation to stem-cell differ- 
entiation and tissue development Jg]. Force is generated in 
cells by myosin motors that pull on the cross-linked actin fil- 
aments comprising the cytoskeleton. Most theoretical studies 
model these forces via the activity of force-dipoles that are 
embededded in a material described by linear elasticity [7]. 
These forces can act within the cytoskeleton resulting in or- 
dering of its actin filaments On a larger scale, the entire, 
contractile cell [9] can be represented as a force dipole that 
deforms its extracellular environment to produce strains and 
stresses that result in effective elastic interactions with other 
cells H lQfl . However, both the cytoskeleton and the extracellu- 
lar matrix comprise cross-linked, semi-flexible polymeric fil- 
aments that respond linearly (with elastic constants that are 
stress independent) only for small stresses. For stresses be- 
yond a critical value, these gels show a nonlinear response 
with power-law stiffening of the elastic moduli with increas- 
ing stress d 1 in . Indeed, nonlinear behavior was measured in 
the elastic and viscous response of biological cells 1121 . 

We report on the first step in theoretical understanding of 
how the nonlinear mechanical properties of an elastic medium 
can radically change the strains and stresses generated by in- 
ternal forces in a long-range manner. Our results are impor- 
tant for the understanding of deformations induced by inclu- 
sions in nonlinear elastic media, active forces in in-vitro acto- 
myosin gels; they also have implications for the interactions 
of acto-myosin segments in cells as well as interactions of 
cells in nonlinear biopolymer media 111311 . We consider an in- 
finite, homogeneous, and isotropic compressible elastic mate- 
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FIG. 1: a) Spherical force dipole generates radial displacement uo at 
radius Rq. b-c) Constitutive relations for simple shear: Eq. (TJ with 
n — —1,6= 100, fi = 1: b) Shear stress vs strain diverges at 7 ma x. 
c) Differential modulus vs shear stress asymptotes to G oc a 13 . 



rial whose response is linear for small stress, but shows strain 
stiffening, with a power-law increase of its differential shear 
modulus with increasing stress, see Fig.QJ. We analyze the re- 
sponse of this medium to a spherical force dipole which exerts 
an isotropic radial force on its surface, see Fig. QJ. This can 
represent a spherical inclusion in a gel and is also motivated 
by the forces exerted by relatively symmetric IU411 contrac- 
tile, non-motile cells such as endothelial cells for which ad- 



hesion sites are uniformly distributed on their perimeter I15H . 
Our model predicts analytic scaling laws as well as numerical 
solutions for the strain that the force dipole generates in the 
medium, as well as its associated energy cost. Our results sug- 
gest that the long-range interactions of such force dipoles is 
significantly modified by the nonlinearity of the medium [16]. 

For simplicity, we first describe the force dipole by the ra- 
dial displacement u that it generates at a distance R from its 
center. We will later invert the problem to obtain the force ap- 
plied at R = Rq. The deformation decays with the distance R 
from the center of the dipole. At R ^> Rq the stress is small 
enough so that the distance dependence of the displacement 
is identical to a linear medium, u(R) m u c s(Rq/R) 2 , but 
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with a coefficient u e g that can differ significantly from that 
of a linear elastic medium for which u e g = Urj. We find that 
for weak nonlinearity (small u ), the deviation from the linear 
solution scales as u e g — u Q oc Uq, whereas for strong nonlin- 
earity (large uq) u e s scales as u e g oc Uq, with the irrational 
exponent 9 = — ^7= w 5.32. These scaling laws induce 
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related, nontrivial scaling relations between the force dipole 
moment, the induced strain and stress, and the total elastic en- 
ergy required to generate the deformation, which are the focus 
of recent traction force microscopy experiments 11411 . 

A gel that stiffens as shear stress is increased may be mod- 
eled by the following elastic energy density functional 1 17]: 



W 



li- 



1 



Hi ~ 3 



- 1 



(1) 




This applies to a system in which the strain is not necessar- 
ily small so that J = det(F) measures the compression with 



Fij = -gjj- the deformation gradient tensor, X the reference 
position and x the deformed position. The shear deformation 
is given by I\ = ii/J 2 / 3 with 1\ = tr(B) and Bij = Fn~Fjk 
the left Cauchy-Green strain tensor 01811 . For vanishingly 
small values of the dimensionless parameter b characterizing 
the nonlinearity, (fl3 yields a compressible neo-Hookean form 
W = % (h - 3) + £ ( J - 1 " ' 

j Fik -1^- oDtains tne rorm <x = 
1)1, with 1 the unit tensor. For small deformations 

Wt + mfi the linear 
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h -y (J — 1), and the Cauchy stress tensor 
obtains the form <x = f B ' 
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1 the unit tensor. 

sa 5ij + 2eij, with e# = \ (j 

strain tensor, and u = x — X the displacement. In this limit, 
I\ pa 3 + 2tr(e) and J « 1 + tr(e), thus the constitutive re- 
lations reduce to Hooke's law, cr — 2/ie + (K — tr(e)l, 
with shear modulus /1 and bulk modulus K. 

For simple shear (x = X + jZ, y = Y, z = Z), (Q~|i gives 
u I2 = fj/y (l + ^7 2 )™ ■ In the limit of small shear, the 
differential shear modulus G = i s constant and equal to 
fi. For n < and b > 0, cr^ diverges as the shear increases 
and the strain 7 approaches 7 max = y/~r> while the shear 
modulus stiffens as G oc with /3 = see Fig.[T]3-c. 

Thus b determines the maximal strain 7 max , which is exper- 
imentally known to be related to the concentration of actin, 
collagen or crosslinkers in biological gels [11]. Prelimi- 
nary estimates based on the experimental data yield values 
2 < b < 400. Theoretically, strong nonlinearity may be intro- 
duced at extremely small strain by taking arbitrarily large val- 
ues of b, which induce arbitrarily small values of 7 max Il9tl . 
Since many biopolymer gels are nearly incompressible, we 
assume the compression is small and treat it linearly, n de- 
termines the exponent j3 that quantifies the strain-stiffening 
behavior of G vs a. The nonlinearity of semi flexible chains 
implies that /3 = | 111 111 , which is obtained by taking n = — 1. 

For the spherically symmetric geometry of Fig. QJ, the de- 
formation is described by the radial displacement u(R). We 
focus on extremely nonlinear materials with 6^1, for which 
the nonlinearity enters at extremely small strain. Here, the 
strain scales as u/R, therefore u/R <C 1 and it is tempting to 



FIG. 2: (Color online) Normalized effective far-field displacement vs 
dimensionless parameter A = bu^ / 7?q quantifying the nonlinearity 
in the system. Color lines indicate numerical solutions of Eq. ([2} for 
different Poisson ratios, and n = — 1. Dotted lines are theoretical 
solutions: Eq. {5]( at small A and Eq. {9]l at large A. 



keep only the terms linear in u. However, since bu 2 is not nec- 
essarily small, we keep both the linear terms and terms of or- 
der bu 3 . We derive the stress tensor from the derivative of (Q]) 
with respect to the deformation, employ this small-strain ap- 
proximation, and eventually obtain the following equation of 
mechanical equilibrium in which the divergence of the stress 
vanishes for R > Rq where there are no internal forces [18]: 
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? \m~ r)> ° 2 = 
h - 3) L - Cx, u = u/uq, R = R/Ro, and A = bul/R%. 
Weak nonlinearity, A <C 1: The linear (A — > 0) solution is 
u(R) = l/R 2 . For Am 2 <C 1, linearization leads to 
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Poisson ratio. For any A this is valid for large enough R. For 
A <tc 1 it is valid already from R = 1, where u(l) = 1. Thus 
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implying that u e s — uq oc Uq. Agreement with the numerical 
solution of Eq. ©maybe seen in Fig.|2]to hold up to A ps 0.1. 

Strong nonlinearity, A 3> 1: The deformation has a shear 
component, I-y — 3, and a compressive component, J. In our 
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model, only the resistance to shear stiffens with increasing 
stress. Hence, for strong nonlinearity, the deformation energy 
will be dominated by shear, so that energy minimization im- 
plies that Ix — 3 should be minimal. A zero shear deformation 
is derived from (0 with 



u = R° 
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and a = ~ 3 + v/ ^ « 0.44. The reason for this irrational expo- 
nent (and for subsequent irrational scaling exponents) arises 
from the homogeneous nature of (0 that results in a quadratic 
equation for a. However, (01 is not an exact solution of Eq. (0 
for all values of R although Fig.0shows that indeed as A in- 
creases, the deformation approaches (0 when R is close to 
1 and the displacement is large. Eq. (ffj) can be satisfied to 
higher order in R by adding a correction to u(R) of order -4 

that causes 1 + ^ ( Jj — 3^ to vanish to leading order: 
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This provides a good approximation for u(R) and in partic- 
ular reproduces the nonmonotonic behavior in the numerical 
solutions, see Fig. [3] In a linear medium, displacements decay 
monotonically. Namely, for a contractile force (uo < 0), the 
radial strain satisfies ^ < for all R > Rq, and ^ > only 
for R < Rq. In a nonlinear medium, 4|| > in a large region 
near the inclusion (R < R*) and ^ < only for large values 
of R > R* where Eq. (01 is valid. One way to understand 
this is that for nonlinear elasticity there is a strong penalty on 
shear strains above a critical level, and extending the region 
in which ^ > helps to reduce the shear strain. Finally, i?* 
represents an effective size of the force dipole, whose strain 
field extends over a distance much larger than its size Rq. 
Defining R* as the position of the maximum of (01, 
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we find that for 4> 1, higher order corrections to u(R) are 
functions of the scaled distance R/R* l20Tl . The increasingly 
good agreement of these corrections with the numerical solu- 
tions as the nonlinearity A increases is seen in Fig.0 For large 
A, i?* ^> 1 or i?* ;§> Rq. The increase in i?* with A suggests 
that even at relatively large distances from Rq, the inclusion 
cannot be treated as a point dipole. This may have important 
consequences for the interactions between force dipoles in a 
nonlinear medium. Our results suggest that these interactions 
will be significantly and qualitatively increased for large val- 
ues of A. Such interactions at surprisingly long distances are 
consistent with recent experimental measurements of biologi- 
cal cells on nonlinear elastic substrates 11311 . 

At long distances the displacements obey u(R) = u e s/R 2 . 
The crudest way to obtain u c ff is by matching this asymptotic 
large-i? form with (0) derived for small R. Thus 
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FIG. 3: (Color online) Normalized displacement vs normalized ra- 
dius for n = — 1, Poisson ratio v = 0.4, and various values of the 
dimensionless parameter A = buy/R^ characterizing the nonlin- 
earity in the system. Colored solid lines are results of numerically 
solving the equation of mechanical equilibrium {2]l. Stars denote the 
locations of the maxima. Dotted blue line is Eq. l|6}, dashed magenta 
line is Eq. 10 and solid black line is higher order correction 1 20] . 



As seen in Fig. [2] the scaling with A is in excellent agreement 
with the numerics while the prefactor differs slightly. Higher 
order terms both in the small- R solution (01 or the higher- 
order corrections given in [20], and in the large-i? solution (0 
yield a more precise prediction for u c g . However, the scaling 
of w e ff with i?* (or with A) remains unaltered, and even with 
such more complicated matching schemes there is still some 
discrepancy in the prefactors compared to the numerical so- 
lution. Noting that -u c ff = u c s/uq, and A oc Uq, we rewrite 
(0 as u cS oc u® with 6 = ^- = « 5.32. It will 

be interesting to test this relation experimentally by inferring 
the near- and far-field displacements proportional to uo and 
Moff from the motion of fluorescent beads in nonlinear media 
deformed by contractile cells. 

Stress: So far we considered the boundary-value problem 
and solved the displacement field u(R) given a certain dis- 
placement uo at radius Rq. The physically and biologically 
interesting questions have to do with the relations between 
the total force AitRqCTo applied at i?o and the deformations it 
generates, where cto = a rr (Ro). In this context, two separate 
questions should be asked. One with respect to the deforma- 
tion uo at Ro and the second with respect to the deformation 
u e ff felt at R > Rq. 

In the linear case (A = 0), cr in = - jj^ oc u = u cS . For 
weak nonlinearity (A <C 1), we find Il20fl that the correction 
to the stress compared with the linear solution is linear in A: 
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For strong nonlinearity (A 3> 1), we find Il20ll that the stress 
increases with A with an irrational exponent related to a: 
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In ( [Tol l and ( fTTT i we used the previously established connec- 
tions between u c s and uq in these two limits. To summa- 
rize, we find nontrivial scaling relations between the applied 
stress cjq and the displacement uq generated where the force 
is applied (R = Rq), as well as nontrivial scaling relations 
between uq and the effective displacement u c ff felt far away 
from the force dipole. However, taken together, these rela- 
tions reduce to a linear relation between oo an d u cS- Thus 
the renormalization of u c g of the linear solution for R 3> R* 
also characterizes the renormalization of the stress applied at 
R = i?o- Interestingly, the same relation cto oc u e s is true 
for both weak (A <C 1) and for strong (A 3> 1) nonlinear- 
ities. Such renormalization is similar to that investigated for 
charged colloids 112 111 

These results may help resolve the question of whether cells 
sense stress or strain [22] by measurements on identical cells 
placed on nonlinear gels with different values of // (and b). 
Since <tq oc /wt e ff, if the cells exert constant stress, one would 
measure u c g oc If, on the other hand, cells exert constant 

strain, Eq. (|9]l leads to u e s oc b 2 ~ 2a . Experimentally, /i and b 
are not independent, but are roughly related by /i oc 7~^ x ^ 



b 2 , with 3 < a < 6 111 1H . This leads to u e g oc fj/i with 
0.7 < T) = ,-, 2+ °L < 1-4, which is easily distinguishable 
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from u e s oc — . 

We now connect our results to additional experimentally 
measurable quantities ffl . and briefly discuss how the total 
energy stored in the elastic deformation (the self energy of the 
force dipole), scales with the force and displacement. In our 
formalism this energy is given by U — J Wd 3 r. However this 
is also equal to the work performed by the active forces start- 
ing from an undeformed state. Since the external forces are 
applied only at R , this simplifies to U = 2ttRq(JoUo oc uoUq. 
In the linear case, U oc a 2 (or Uq). For strong nonlinear- 
ity ( A > 1), we obtain U oc cr^ +1/e (or u Q +1 ) w afr 19 (or 
Uq 32 ). A direct comparison with the nonlinear scaling laws 
deduced from traction-force microscopy experiments 11411 is 
beyond the scope of this paper, since these measured displace- 
ments assumed linear elasticity for calculating force and self 
energy. Moreover, in the above analysis we assumed a fixed 
spatial size Rq for the force dipoles and studied the depen- 
dence only on the dipole strength as parameterized by uq or 
<7o, whereas in a more realistic model there may be correla- 
tions between the radius of the spherical dipole (or cell size) 
and its activity. Nonetheless, our predicted irrational scaling 
exponents should show up in such measurements. 

In summary, we predict scaling laws for the response of a 
nonlinear elastic medium to an internal, spherical force dipole. 
We find that the nonlinearity changes the induced strain from 
one that is constant in sign to a strain that changes sign at a 
distance i?„ that for strong nonlinearities can be very large 
compared to the size of the dipole. We also predict scaling 
laws for the displacement that behaves as a power law with 
irrational exponent. Moreover, we identify nontrivial scaling 
relations between the applied force (do) and the deformations 
generated at the dipole surface (mq), and between these de- 



formations and the deformations felt at long distances (u e g). 
However, surprisingly, we find that the stress applied by the 
spherical dipole can always be written as a linear function of 
the renormalized displacement u c g characterizing the strains 
induced at very large distances. The magnitude of the ef- 
fects we find can significantly enhance the interactions be- 
tween force dipoles. This will have important implications 
for cell-cell interactions in nonlinear elastic media as already 
suggested by the experiments in II 1 3 1 - 
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